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We present analytical calculation of the first seven moments of the heavy quark 
vacuum polarization function in the three-loop order. The obtained results are 
compared against the asymptotic formulas following from the threshold singulari- 
ties. We also discuss the fi dependence of the moments within the BLM procedure. 



1 Introduction 

The quark vacuum polarization function 

Uf{q2) = 3?/ dxei9X{ ° 1 T tiWrtW |0) ' (1) 

with jf = ifi fipiftf being the vector current for a quark /, is an interesting 
object related to a number of important physical quantities. An (incomplete) 
list includes: 

• The combination^] 

D ® = ? T+d™m with nem = E<ffn' (2) 

is the quark contribution (to order e 2 ) to the photon propagator D(q 2 ). 

a We ignore the so-called singlet contribution proportional to (^2f Q/) 2 ant ^ ^ rs * spear- 
ing in order orj. 
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• The optical theorem relates the inclusive cross-section of e + e~ annihila- 
tion into hadrons and thus the function R(s) = a tot I '° 'point to the discon- 
tinuity of n em in the complex plane 

R(s) = 127rImir m (s-He). (3) 

Conversely, the vacuum polarization is obtained through a dispersion 
relation from its absorptive part, vis., 



s(s — q 2 ) 

• Replacing the electromagnetic charge Qf in H em by the vector weak 
coupling Vf = 2/3 — 4Q / sin 9 W one arrives at the vector part of the 
quark contribution to the Z boson propagator and, again through the 
dispersion relation, to the vector part of the decay rate of the Z boson 
to hadrons. 

To order a s the calculation of IL£ was performed by Kallen and Sabry in 
the context of QED a long time ago a. With measurements of ever increasing 
precision, predictions in order a 2 are needed for a reliable comparison between 
theory and experiment. _ 

Such a calculation was recently presented inl3. It employed a semianalytic 
approach based on conformal mapping and Pade approximation I3qEIi3. Im- 
portant ingredients were the leading terms of the expansion of IF in the high 
energy region (—q 2 )/m 2 ^> iQ'EI, the Taylor series around q 2 — which has 
been evaluated up to terms of order (q 2 ) and information about the structure 
of II in the threshold region. The method was tested in the case of the "double 
bubble" diagrams against the known analytical resultEl with highly satisfactory 
results. The calculation leads to a point- wise prediction of the absorptive part 
W (s) in order a 2 with conservatively estimated accuracy of about 6%, which 
should be sufficient for comparison with experimental results in the foreseeable 
future. 

An important case when the full mass dependence of IF is necessary is 
a precise determination of a s and m from QCD sum rules for bb E£l. This is 
because of the still persisting contradiction between the values of a s determined 
from low- and high-energy measurements c3. Fortunately enough, QCD sum 
rules for bb are based on the moments of production cross section of bb states 
in e + e~ annihilation which can been computed exactly even in order a 2 . 

In the present paper we report on analytical calculation of the first seven 
moments of the function or, equivalently, of its Taylor series around q 2 = 
up to (and including) terms of order (q 2 ) 7 . 
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2 Notations 



We deal with the case of QCD containing n; = Nf — 1 massless quarks and 
a massive quark ipp with the (pole) mass mj? later referred to as m. The 
corresponding polarization function Tl F and the ratio R F will be denoted as 
IT and the ratio R, respectively. It is convenient to define 

ii( q 2 ) = n<°' + ^c f nw (5) 

n (2) 



c 2 F nf + CaCfU^I + c F T ni n[ 2) + c F rn {2) 



with a s {jx) being the MS coupling constant. For every polarization function 

with i = 0, 1, 2 and ? = A, NA, I, F (if i = 2 ) we define its moments C ? (i) 
as follows 

n>0 v ' 

or, equivalcntly, 

with t> = y^l — ^f-. Note that by definition 

= C 2 F C {2 1 + C A C F Cf Ajn + C P TmC$ + C F TC i2 l (8) 

In order to fix the absolute normalization we put below the lowest order result 
for R(s): 

R(°)(s)=3 1+ U/l . (9) 



3 The calculation 

Important information is contained in the Taylor series of n(g 2 ) around zero. 
The first seven coefficients of the Taylor series around q 2 — are calculated 
with the help of the program MATAD which has been used previously for the 
three-loop corrections to the p parameter!! 2 ], to ArEJ and the top loop induced 
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reaction e + e — > Hl + l -m. The package is written in FORMEJ. For a given 
diagram it calculates the derivatives with respect to the external momentum 
q 2 to the desired order, raising internal propagators to powers up to about 
20. It then performs the traces and uses recurrence relations based on the 
integration by parts method t£EZI to reduce the resulting three-loop tadpole 
diagrams (in dimensional regularisation) . In order to save space we present 
our results converted in numerical formJn the Table 1. The original analytical 
expressions can be found in our paperE3. 

The coefficients c} 2) - cf ] were already presented ini, cj?\cj^ and 
Cy^-are new. The C\ term of — cjp and are in agreement with! 
andEfl, respectively. The results in the Table are given for a particular choice 
of \l = to; in order to transform them to a general /i the following relation 
should be used: 

CL 2) (//) = C< 2 > (/*) + ClPp ^ (-) 2 with A) = Y2°A - jNj (10) 



n 


r {0) 




r (2) 
U A,n 


r (2) 

NA,n 




r (2) 


r.(2).s 




1 


1.067 


4.049 


5.075 


7.098 


-2.339 


0.7270 


31.66 


33.22 


2 


0.4571 


2.661 


6.393 


6.311 


-2.174 


0.2671 


30.99 


32.44 


3 


0.2709 


2.015 


6.689 


5.398 


-1.896 


0.1499 


28.53 


29.79 


4 


0.1847 


1.63 


6.685 


4.699 


-1.671 


0.0995 


26.29 


27.40 


5 


0.1364 


1.372 


6.574 


4.165 


-1.494 


0.0723 


24.41 


25.41 


6 


0.1061 


1.186 


6.426 


3.746 


-1.353 


0.0557 


22.84 


23.74 


7 


0.0856 


1.046 


6.267 


3.409 


-1.239 


0.0446 


21.5 


22.33 



Table 1: The results for the first seven moments. The last two columns display the moments 

(2) 

Cn for the cases of Nf = 5 and JV^ = 4, respectively. 



4 Higher n moments and threshold singularities 

In this section we compare our exact moments with their asymptotic behaviour 

(i) 

originating from the corresponding threshold singularities, with r) n standing 
for the ratio C^JG^ as . 
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As it follows from (]?]) the asymptotic behaviour of the moments in the limit 
of large n is completely determined by the behaviour of R(s) at the threshold 
s m Am 2 . The constraints on the threshold behaviour of Tl(q 2 ) originate from 
our knowledge about the nonrelativistic Greens function in the presence of 
a Coulomb potential and its interplay with "hard" vertex corrections. For a 
theory with nonvanishing (3 function (QED with light fermions or QCD) the 
proper definition of the coupling constant and its running must be taken into 
account. Below we cite results that are available in the literature about the 
threshold behaviour of R(s) and its moments. 

R {0) - thr = a(y-y). (U) 

~ n 3/2 8n6 /a- W 
R^> thr = 3^tt 2 -6«V (13) 



C\ part 001: 
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i 

n n 3 / 2 



C F C A parti: 



8v 

r 9/2 

Qn 1 / 2 n 



R f'*>» = 3 f^-STr 2 ), (15) 



= ^-4-. (16) 



C F riiT parti: 



^l(H_ llln4 + 11 vI N 



= —I- n-' iO.-l- u*> 



i?f>^ = 37 r 2 (±ln^-^), (19) 




Cg 1 " = -|-^+ ln4 -^7]- ( 2 0) 
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CpT parti 



(22- 27r> + (-^ + ^r> 3 , (21) 

In Table 2 we compare the asymptotic formulas with the available mo- 
ments. We observe that at n = 7 the first two terms of the asymptotic expan- 
sion in 1/n agree with the exact results with the accuracy of about 35%. 



n 


r W 

' n 


' n 


'A,n 


1 NA,n 


(2) 
r l,n 


(2) 
r F.n 


-W.6 
' n 


J2),4 
' n 


1 


4.814 


1.457 


-0.474 


-1.755 


2.536 


2.51 


-0.845 


-0.902 


2 


1.297 


1.096 


10.51 


2.522 


1.032 


1.28 


5.522 


4.624 


3 


1.121 


1.046 


1.937 


1.700 


0.9709 


1.128 


2.078 


1.982 


4 


1.067 


1.031 


1.479 


1.499 


0.9611 


1.075 


1.643 


1.597 


5 


1.043 


1.024 


1.322 


1.407 


0.9608 


1.049 


1.472 


1.442 


6 


1.03 


1.020 


1.243 


1.353 


0.9628 


1.035 


1.38 


1.357 


7 


1.022 


1.018 


1.197 


1.317 


0.9653 


1.026 


1.321 


1.304 



Table 2: The results for the first seven ratios of exact moments to their asymptotic values. 
The last two column display the ratios for the case of Nf = 5 and Nf = 4, respectively. 



R 



(2),thr 



C 



(2), os 



i0| ) is to the running of a s (/it 2 ). 
to our results. It suggests to 



5 Choice of fi 

The /j, dependence of C^- 2 ' as displayed by eq 
Now we want to apply the BLM procedure 
choose the scale of the 0(a s ) term in such a way that the contribution of the 
0(a 2 ) part proportional to f3o is absorbed. This prescription is based on the 
observation that the remaining coefficients of the a 2 s terms are often relatively 
small. It is possible to treat each term of 0{z n ) separately. In Table || we list 
the BLM scale [iblm for each Taylor coefficient together with the numerical 
value of the original (C^ 2 ') and the 0{a 2 s ) term which remains after the BLM 
scale is adjusted (C^ LM ). It is interesting to note that \xblm is decreasing 
with increasing n. This is plausible because the higher Taylor coefficients are 
increasingly dominated by the threshold region and the characteristic scale 
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n 


r {0) 


r (i) 


r (2) 






rW 

On 


My 

On 


1 


1.07 


4.05 


30.1 


13.7 


0.420 


1.92 


3.82 


2 


0.46 


2.66 


29.5 


14.3 


0.294 


0.83 


3.69 


3 


0.27 


2.01 


27.3 


14.0 


0.244 


0.39 


2.50 


4 


0.18 


1.63 


25.2 


13.5 


0.215 


0.15 


1.65 


5 


0.14 


1.37 


23.4 


13.0 


0.195 


0.01 


1.12 


6 


0.11 


1.19 


21.9 


12.5 


0.181 


-0.09 


0.85 


7 


0.09 


1.05 


20.7 


12.0 


0.169 


-0.15 


0.76 



Table 3: The BLM scale fiBLM is expressed in terms of the original scale m. For the 

(2) "(2) _ 

numerical values of C„ and C 7l n; = 5 and /i = m, respectively, (i = m is used. The 

double bubble diagrams with two massive fermion loops are also included. C^ LM is the 

coefficient remaining after the /3q term is absorbed. 



at threshold is the relative three-momentum of the quarks. Note that C„ 

(2) 

remains nearly constant whereas C n decreases for increasing n (but remember 
the rapidly decreasing coefficients Cn" 1 )- 

In comparison to C n 2 ^ in Table |^ also the corresponding one- and two- loop 
coefficients (C n °\ C'rP) and the two- and three- loop coefficients (Cn \ C n 2 ^) in 
the MS scheme are listed. Whereas Cn^ is roug hly of 0(1) C n 1] varies from 
approximately 1.9 down to 0.01. Also the sign changes. Therefore the BLM 
procedure is rather unstable and not applicable. On the other hand, the MS 

~ (2) 

coefficients C n are already reasonably small, so there is no urgent need for 
an optimization procedure. 



6 Conclusions and summary 



We have presented analytical results for first seven moments of the heavy quark 
vacuum polarization function to order 0{a 2 s ). We found that BLM procedure 
does not lead to a significant improvement of the perturbation theory for the 
moments. It is demonstrated that the threshold singularities do describe the 
higher moments (with n = 6, 7) with inaccuracy not exceeding 40% (in order 
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